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ABSTRACT 

We present stationary models of protostellar population iii (Pop iii, for short) accretion discs, 
compare them to Pop i discs, and investigate the influence of the different chemical com- 
positions on the occurence of gravitational, thermal and thermal-viscous instabilities in the 
discs. In particular in the cooler regions, we find major differences between Pop iii and Pop i 
discs, both in the structure and stability behaviour This is mainly due to the absence of most 
molecules and dust in Pop iii, which are very efficient absorbers in Pop i discs. 
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1 INTRODUCTION 

The important role of angular momentum during stellar and plan- 
etary formation is best evidenced by the remarkable co-planarity 
• of the planetary orbits in our solar system. Observatio ns of proto- 
stellar and protoplanetary accretion discs (Hartmann 1998) point 
towards this being generally the case. Even in the lack of direct ev- 
idence it is reasonable to assume that — as far as the role of angular 
momentum is concerned — the situation was not substantially dif- 
ferent for the very first generation of stars, the population iii (Pop 
III, for short). 

While the general role of angular momentum may have been 
the same, and certainly will have led to the formation of popula- 
tion III accretion discs, the structure and the stability of these discs 
may have differed considerably from that of contemporary proto- 
stellar accretion discs. This difference in structure and evolution is 
not surprising given the difference in the chemical composition of 
accretion discs now (Pop i) and then (Pop iii), in particular the lack 
of efficient cooling agents in the primordial matter discs. 

Models of accretion in early universe so far assumed spherical 
accretion I Stabler et al. 1986). Only recently pe ople starte d to wor k 
on accretion through a disc (Tan & McKee 2004, herafter lTM04b . 
In particular, the importance of the opacity coefficient for th e stabil- 
ity of accretion discs h as been known for quite a while (e.g. lPuschll 
ll993llHure et al 1200 1*). In this contribution, we wish to investigate 
the structure and stability of population iii accretion discs. For this 
purpose, we present stationary, vertically integrated, axially sym- 
metric discs models. 

In Sect.|2|we shall summarize the properties of the accretion 
disc model we used and present the set of equations solved. The 
subsequent Sect.|^is dedicated to a short summary of the numerical 
method used to solve these set of equations. We then (Section |4} 
discuss the opacity coefficient for a Pop iii chemical composition. 
Before we proceed to the presentation of the accretion disc models 
(Sect. 0, we shortly review criteria for thermal, thermal-viscous 



and gravitational instabilities in Section |5| and give formulae for 
the timescales present in discs (Sect.|6} developing handy formulae 
for relevant limiting cases. We end in Sect. |H| with discussions of 
these results and put them into perspective. 



2 THE ACCRETION DISC MODEL 

We use a standard one zone accretion disc model B Pringlelll98ll) 
with slight but i mportant modifications compared to the results 
of iMaven i2003h : Our discs are allowed to be vertically self- 
gravitating and can become optically thin. Vertical selfgravity 
means that in a direction perpendicular to the rotational plane of 
the disc the gravitational acceleration is no longer restricted to be 
due to the projected contribution from the central object only, but 
now can also be due to the local mass distribution, or both at the 
same distance from the central object. The gravitational pull in ra- 
dial direction is still only due to the central force as we exclude so- 
lutions where the disc mass becomes larger than the central mass, 
i.e. the rotation al law is still Kepler ian (KSG discs, Keplerian self- 
gravitating, see IPuschl et alj200(]l) . 
We have 

• the continuity equation 

M = InrlLur 

• the angular momentum equation (v viscosity; f{r) accounts 
for the inner boundary condition, see eq.0 
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the energy equation 
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It has been shown JTMOJ) that the ionisation energy may not be 
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neglected although the thermal energy is negligible. We do not ac- 
count for this effect in our models (see discussion in Appendix IdI. 
• hydrostatic equilibrium in vertical direction 
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The pressure due to matter and radiation is sustained by the vertical 
component of the gravitational force g;. h of the central body and the 
vertical selfgravity of the disc (^ksg = 2;rGX, Paczvnski ( 1978)) 

The optically thick -^ optically thin transition is achieved 
by modifying the e quation of state and the optical depth 
JArtemovaetalll99el) 



kT 4cr ., 
P = Pg+Pr = p + r-* tr + 

The effective optical depth is given by 
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where we define the Rosseland and Planck optical depth by 

Tr/p = -2*-R/p(p, T) 



(4) 



(5) 



(6) 



The radiation pressure reaches the black-body value p = ^j'* 
for Tgjj = tr — > oo and is negligible in the optically thin regime 

4 

(tr < 1, Tgjj = -^ — > oo, 2.^ —, 0), In a optically thin medium 

photons are not able to contribute to the pressure in the gas as the 
mean-free path becomes long compared to the physical size of the 
medium they are passing through. The formal optical depth in eq. 
J5j which enters into the energy eq. J2j equals the usual diffusion 
approximation for tr » 1 and moves into the optical thin regime 
for Tr <c 1 . In this respect we account for the transition from the 
optically thick diffusion approximation (Q^ = ^T*) to the opti- 
cally thin emission {Q^ = 2crT'*Tp. 
The inner boundary condition (Eir 



: 0) can be expressed by 
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Once we are sufficiently away from r,, f{r) 
we set /(r) = 1. (see Sect. lT.l.n 

We furthermore define 2 parameters 

Pp - — — T}VS,G 



gc.H 



(7) 

1. For our models 

M 



Pg + Pr 



gc.H + gKSG M + Anpr^ 



f}p measures the contribution of the gas pressure to the total pres- 
sure (ySp = 1 : pure gas pressure, j8p = 0: pure radiation pressure), 
and ;/KSG the vertical selfgravity (77KSG = 1 : no vertical selfgravity, 
jyKSG = 0: fully vertically self gravitating). 

For the viscosity prescription we use a-viscosity 
JShakura fesiinvaev 1973) 

y = a^ (8) 

This corresponds to a stress tensor element t,^ = -UpP with ap = 
I a which translates into the more commonly used v = acsH for 
the non-selfgravitationg case. We use this more general viscosity 
prescription to avoid the difficulties of isothermality and rapidly 
increasing density in the vertically and fully selfgra vitating parts 
of accretion discs when using v = acsH in all cases iPuschl et alj 
L2OOO) '. 

' We note that in lPuschl et alji200(ll 's notation the viscosity parametriza- 



M 


accretion rate 




E 


Surface density 1 = 2pH 




P 


density 




T 


temperature 




H 


disc half thickness 




V 


viscosity prescription 




It, 


Radial drift velocity 




n 


Rotation frequency Hnsg/ksg = 


_ IGM 


p 


pressure 




Cs 


sound speed 




^eff 


effective optical depth 




M 


central mass 




r 


radial distance 




*fR/P 


Rosseland and Planck opacity 




P 


mean molecular weight 





Table 1. Nomenclature of the variables used. If not stated otherwise, the 
symbols refer to the disc's mid-plane 



3 NUMEIUCAL TREATMENT 

As we do not treat radial selfgravity and therefore have no implicit 
radial coupling, we are able to reduce the equations above into two 
equations depending only on density and temperature. 



\4pr3 j\^nv 



kT AcrT^ I KrM 4 



3 GMM 



jimp 
Act 



3cTeff \ 67ry 3 



Sfeff 






(9) 



(10) 



We solve these equations using a 2D nested intervals method 
for (p, r) in the range of -16 < logp < -2 and 1.8 < log T < 4.6. 
We refine the initial rectangle 10 times so that we end up with an 
accuracy of A (logp) w 0.02 and A (log T) » 0.003. 

In our system of equations we have left the effective optical 
depth Tgff (given by eq.|5j and the viscosity v. For the cc- viscosity 
prescriptions it is impossible to explicitly calculate y (p, T) as the 
radiation pressure depends on the surface density which itself de- 
pends on the viscosity. This implies a polynomial equation for v of 
order 3. The details of this calculation are given in Appendix IXI 
The optical depth still containing the surface density can be fully 
expressed in terms of density and temperature by using the angular 
momentum equation Q. 



4 OPACITY 

Primordial matter as an unprocessed relict from the big bang 
strongly differs from today's chemical composition as then the 
heavier elements which have to be produced first in stars were not 
present. Due to the lack of these metals, no dust could form. This 
lack of dust affects the opacities at low temperature, as dust is a 
very strong absorber. 

In the absence of good absorbers unusual absorption effects 
become important: in the high-density and low-temperature regime 
collision-induced absorption (CIA) of H2/H2, H2/He, H2/H and 
H/He-pairs dominates. In the low-temperature and low-density 
regime we have Rayleigh scattering in the Rosseland opacity while 
CIA continues to be the dominating effect in the Planck mean. 

As there was no consistent and complete set of Pop iii opaci- 
ties available in the literature (except the ones of Lenzuni et al. from 

tion J8) corresponds to what they call the dissipation-limited case, i.e., 
where the turbulent velocity is limited by the sound velocity. 
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Pg » Pr 


TR» 1 


Ar - (1 + r;KSG)BR + 3(1 + ;;ksg) < 


Tp <K 1 


-Ap + (1 + )?ksg)Bp + 3(1 + ;7ksg) < 


Pr » Pg 


TR» 1 


4Ar - ;;ksgBr - 4 < 


Tp <K 1 


-4Ap + !;ksgBp - 4 < 



Table 2. Conditions for thermal instability 



Pe :» Pr 



Pr » P? 



TR :» 1 



Tp « 1 



TR :» 1 



Tp « 1 



3Ar -( 1 +fjKSG)fiR+5( 1 +';kSg) 
-*R -( 1 +';ksg)Sr+3( 1 +i;ksg ) 



<0 



-3Ap+(l+;;KSG)fip+3(l+'jKSG) 

-Ap+(1+>;ksg)Bp+3(1+';ksg) 



<0 



12AR+(2-i;KSG)fiR-4(l-2>;KSG) ^n 
4Ar-i;ksgBr-4 



-12Ap+(2;,KSG-l)gp-4(2->;KSG) 
-4Ap+i;ksgSp-4 



<0 



1000-7000 K and OPAL and OP ones above 10 000 K), we com- 
puted them on our own (Mayer, Diploma thesis; Mayer & Duschl, 
in prep.). The results are Rosseland and Planck mean opacities for 
-16 < logp < -2 and 1.8 < logT < 4.6 assuming chemical equi- 
librium and only including continuum absorption (no lines). Our 
opacities match those of Lenzuni better than 10 % above 3000 K. 
Below 3000 K deviations are due to the implementation of CIA 
(e.g. Borvsow et al. 1989)^. The consistency with the OPAL and 
OP opacities is of the same level of accuracy, although we see larger 
differences in the high-density limit due to neglected pressure ion- 
ization. 

During the last years it became widely accepted that it is suffi- 
cient to take hydrogen and helium species into account when deal- 
ing with Pop III matter. However there is a small number fraction 
of deuterium and lithium present, the latter of which must not be 
neglected. 

Lithium belongs to the group of Alkali metals and can there- 
fore easily be ionized. We included a simple deuterium and lithium 
chemistry. The deuterium chemistry does not affect the opacity at 
all, but the lithium changes the opacity up to two orders of mag- 
nitude due to its ability to provide free electrons at comparatively 
low temperatures. We do not include absorption effects of lithium 
and deuterium explicitly, lithium is only needed to alter the chemi- 
cal equilibrium and therefore the absorption properties. The largest 
positive differences relative to the Z = calculation occur at tem- 
peratures around a few 1000 K and at high and low densities while 
there is a negative difference at intermediate densities due to the 
destruction of Hj. 

Our opacity calculations show a difference between true Pop 
III and zero metallicity matter in this temperature range. Currently 
we only included continuum absorption of hydrogen and helium 
species with the exception of Hj bound-bound transitions. Further 
details, a discussion of the opacity calculation and the inclusion of 
line absorption can be found in Mayer & Duschl (2004, in prep.). 
The models presented here were calculated with a Lithium fraction 
of /u = 5 • 10"'" and a Deuterium fraction of /n, = 4 ■ lO"*". 

To compare with Pop i opacity we use the mos t recent compi- 
lation of low-temperature dust and gas opacities bv lSemenov et alJ 
J2003h. 



Table 3. Conditions for viscous instability 
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Ice grains 


(I+J/ksg) 
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Evaporation of ice gr. 


10(1 + ;;ksg) 
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Metal grains 


|(1+'/ksg) 
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Evaporation of metal gr. 


1 + 27(1 + riKSG) 


3+29(1+,,ksg) 
I+27(1+;7ksg) 


Molecules 


2 
3 


6 + 3;;ksg 


H-scattering 


3-7(1 +;;ksg) 


3-15(1+wsg) 
1-21(1+,KSG) 


Bound-free and free-free 


1 + 5(1+'7ksg) 


6+5(l+r,KSG) 
2+(1+i;ksg) 


Electron scattering {fip = 0) 


-4 


1 - 2f;KSG 



Table 4. Stability criteria (eqs. I12|and| 13l evalu ated for the fits to Pop i 
Rosseland mean opacities from Jfiell & LiiJl994B . 

prone to axisymmetric instabilities leading to fragmentation. The 
true causes for low-2 discs are still under debate (e.g. see discus- 
sion in .Goodman-2003.) . 



5.2 Thermal instability 

Suppose we got an equilibrium solution from egns.lgland llOl Then 
the disc is thermally unstable if an increase (decrease) of the heat- 
ing rate is answered by a slower increase (decrease) of the cooling 
rate. In mathematical terms this criterion is given by 



dlogT 



dlogT 



>Q 



(12) 



If we approximate the opacity by a power law depending on 
density and temperature k = Cp^T^, we can express the thermal 
instability condition <12> for limiting cases in terms of quantities 
depending on material functions only. We give the exact derivation 
in Appendix|n] We summarize the most useful formulae in Table|5| 
(All cases for/j = const). 



5 STABILITY 

5.1 Gravitational instability 

The gravit ational instabil ity is usually characterized by the Toomre 
parameter j%om3l963): 



Q 



'ttGI. 



(11) 



where we make use of k- = 4fi^ + 2RQ^ = Or being the epicyclic 
freqeuency for Keplerian rotation. If 2 <c 1, then the disc might be 

^ For an up-to-date reference of available CIA data see 
|http : //www ■ astro . ku . dk/f abory sow/programs/ 



5.3 Viscous instability 

A disc is viscously unstable if an increase (decrease) of the local 
accretion rate is answered by a decrease (increase) of the surface 
density. Mathematically we can write this instability condition as 



aiogM 



<0 



(13) 



51ogZ 

We call a disc thermal-viscously unstable if it is thermally and 
viscously unstable. 

For a power-law approximation of the opacity /cr/p = 
Cop'''*/T^K/p and fi = const we get the results shown in Table|^ 

As an examp le we apply our instability criteria to the fits of 
iBell & LiiJ il994h for Pop i chemical composition with tr » 1, 
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Figure 1. Sample accretion disc model for Pop m (left, Mayer & Dusclil, 2004, in prep.) and today's (right, ISemenov et all2003l) chemical composition for 
M = lO~^M0yr~'. Note the different colouring schemes mandated by the vastly different opacity range. 



assuming gas pressure to be dominant {fip = 1) except for Thomson 
scattering (Jip = 0), and show them in Table|4] 

We see a thermal instability due to H-scattering and a 
thermal-viscous instability for electron scattering (for t/ksg > 
1/2 only) This is the so called Lightman-Eardley instability 
JLJghtman & Eardlevll974h . 



6 TIME SCALES 

6.1 Viscous, thermal and hydrostatlcal equilibrium 

The viscous, thermal and hydrostatic timescales are given by 
IPringldilQsih . 
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6.2 Chemical equilibrium 

A cautionary remark when dealing with Pop in matter should be 
made regarding the chemical equilibrium. Due to the inefficient 
cooling compared to dust (The primary coolant being H2), the only 
way to form H? in the low-density regime is via the H~ and H^ 
channel i McDowelJ 196?: 'Pee bles & Dickell96i:ISaslaw & Zinovl 
[1967). However in the high-density regime (p > 10"'* g cm"^), 
three-body reactions rapidl y convert neutral atomic into molecular 
hydrogen JPallaet alJll98.l) . A review des crib ing chemical reac- 
tions in the early Universe can be found in lLep p et al. ( 2002). We 
calculated the H2 formation timescale using the reaction rate ^4 of 



IPalla et alJ il983h (see Appendix O. We get in the disc geometry 

_7 

,10-'*gcm-3J \^/ ll/ 
with^(M,p,.)' = 1+4.76(^)(^^P,^^^)"(^)" 

7 POPULATION III ACCRETION DISCS 

7.1 Overall Properties 

As we have seen in section |4| Pop in matter gives rise to unusual 
absorption mechanisms like CIA. Due to the lack of good absorbers 
it also cannot cool below a few 100 K (for an ove rview of the ini- 
tial c onditions used for collapse calculations see iRipamonti et alJ 
l2002h with the only efficient coolant being mainly H2. This mini- 
mum temperature compare d to the today's value ar ound 10 K im- 
plies larger Jeans masses iClarke & B romm 2003) and therefore 
larger accretion rates. Some authors conclude therefore that Pop m 
stars need to be more massive. 

Fiducial values for primordial accretion discs are taken to be 
M = lO^^'Moyr"' compared to M = lO^'Mgyr"' for today's accre- 
tion discs. 

In our stationary accretion disc model the only parameters we 
can adjust are the central mass, the accretion rate, and the viscosity 
prescription. As standard value for a we use a = 10"^. 



7.1.1 Self-similar solutions 

For all models presented there is no specification of an inner bound- 
ary. In the scope of this paper, we want to keep the discussion as 
general as possible. 
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Equations J9) and <10> are two equations for the determina- 
tions of the density p and the temperature T at a given radial dis- 
tance r, central mass M and accretion rate M. However M and r are 
coupled in the sense that the equations yield the same result pro- 
vided that n^ = ^, the rotation frequency is the same. Thus we 
present our disc models in terms of a generalised radial coordinate 



log q = log 



\1 AUJ' 



1 



■ loe 



This approach only applies for NSG and KSG solutions. For 
the FSG(full-selfgravitating) case O also depends on the mass dis- 
tribution of the disk itself. 

For all models we choose a radial grid extending from 
I 
\0^^ . . . IO'Mq^ AU. At some accretion rates the inferred radius of 

the central s tar can be larger than o ur formal inner disc radius (e.g. 
see Fig.l in lOmukai & Pallal200lh . The main scope of this paper 
is more on the global properties of Pop iii accretion discs in com- 
parison to Pop I discs. Therefore we do not specify any boundary 
conditions. This even makes the disk models applicable to primor- 
dial BH accretion. 

Sample accretion disc models are shown in Fig.^ 

7.1.2 High-temperature limit 

For high temperatures (log T > 3.8) the models are nearly inde- 
pendent of metallicity, because at high temperatures the opacity is 
mainly dominated by hydrogen and helium species (up to an order 
of magnitude). And due to the T'^/r^ff oc MQ.^ proportionality a 
change in r ex /cr by one order of magnitude changes the tempera- 
ture only by a factor of 1.8 provided E remains constant. 

7.1.3 Isothermal outer disc 

One striking feature in the temperature T(r) plot of Fig.|2|is the 
isothermal behaviour of the outer parts of Pop in discs. The tem- 
perature seems to be locked to approximately 3000 K. Inspect- 
ing the energy equation J2} and replacing Tgff by -^ with Tp = 

^^Cr(^r(^r-eget: 



9ak 



\-Ap\ 



Sa-TgCpfimp 



For all absorption effects the exponent Ap remains smaller than 
unity but the temperature dependence is far stronger. This is be- 
cause at temperatures of a few 1000 K, H" absorption is dominant. 
In the optical thin part (where Planck mean values apply) this ab- 
sorption has a temperature dependence 5p = j^y ~ 20 ... 30. 
The high value of Bp ensures that then the exponent of the last equa- 
tion is nearly zero and thus all radial dependencies vanish. The disc 
is isothermal in radial direction. 

Such a temperature lo cking has also been shown for the spher- 
ical collapse calculations of lStahler et alJ ll986^ in the photosphere 
of the accretion shock of an evolving primordial star. Although the 
mechanism (//" absorption) is the same, the isothermality for disk 
accretion is the property of the mid-plane temperature, while for 
spherical accretion this is evident in the photosphere of the evolv- 
ing star. 

Approximating the H" opacity (neglecting density depen- 
dence, Ap = 0) we have 



with Cp = lO-'**^ and Bp = 22. 

Following our argumentation above we get 



\O.Ol/\MoJ \AU/ 



fp.H 



= Ci 



\3000K/ 



T = 3440K ■ 



In today's accretion discs at 3000 K molecular opacities of 
H2O and TiO soften the temperature dependence and therefore let 
the temperature decrease significantly. With decreasing tempera- 
ture we get into the region of dust opacity which does not provide 
these steep temperature gradients. 

7.1.4 Disc thickness 

As one consequence of this radial isothermality in Pop m discs the 
relative disc thickness H/R increases as the disc is optically thin 
and H/R = cj {nr){l + t^ksg/ d -r]Ksa))^ '^ c,/ (flr) ex (Tr)^-. 
For constant T we have H/R ex r? . This leads to the break-down 
of the thin-disc approximation as H/R reaches and exceeds unity. 
The other consequence is the radial decline of the surface density 
1, cc r^i . We use the angular momentum equation Q and get c^E oc 
SIM and with isothermality S oc Q ex r"i . 

In the dust opacity disc models the temperature is not constant 
in the outer disc. Therefore H/R decreases with the onset of vertical 
selfgravity, as then H/R oc T/ (Zr) » T/r. 



7.1.5 Minimal temperature 

Between the hot inner part and the isothermal outer part, stationary 
Pop III discs have an intermediate region, where the disc can cool 
below the 3000 K of the outer parts. This is mainly due to the onset 
of CIA absorption. But as the density decreases, CIA absorption 
gets less and less efficient (absorption coefficient proportional to 
the product of number densities of the contributing absorbers). At 
a certain point the absorption is so inefficient that the disc becomes 
optically thin and the absorption turns into optically thin emission. 
The temperature has to rise again until it reaches that of the tem- 
perature locked outer region. The minimal temperature decreases 
with increasing accretion rate as then more and more material is 
accreted. Therefore the densities get higher and so does absorption. 



7.2 Instabilities 

7.2.1 Gravitational instability 

Figure 121 shows the Toomre parameter for Pop m and Pop i discs. 
While the Toomre parameter decreases monotonically for Pop i 
discs (as temperature does) there is a difference in the Pop m case. 
The temperature reaches a minimum value and so does the Toomre 
parameter. 

With the analytical considerations of Sect. l7.1.4l we can write 
the Toomre parameter <1U for the optically thin isothermal outer 
region as 

\0.0l/\l0-4Moyr-V \3440K/ 

For accretion rates in excess of 1.07 ■ 10^'Moyr"' (a = 0.01) 
the necessary condition for fragmentation is fulfilled in the Pop m 
case. 

In isothermal regions the Toomre parameter depends only on 
radially constant quantities. The surface density E and the rotation 
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Figure 2. Comparison of the properties of discs with primordial and today's chemical composition: mid-plane temperature T{r), surface density S(r), Rosseland optical depth T]i(r), relative thickness H/r, and mean 
molecular weight /j(r). Solid lines are stable solutions, dotted lines coirespond to thermally unstable solutions and dashed lines are vertical selfgravitating solutions. Dash-dotted lines are vertical selfgravitating 
thermally unstable solutions. 
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Figure 3. Toomre parameter Q for a accretion rate of 10 ''Moyr ' in disc models of Pop in (left) and Pop i (right) material. The disc is stable with respect to 
fragmentation for log Q> 0. 



frequency Q. both show the same radial dependence which cancels 
in eq. <11> . 

The maximum value of the viscosity parameter for 
marginally stable, selfgravitating discs is still a matter of debate: 
[Tan & Blackman (2004J use a = 0.3, following Goodman 1 2003') , 
while Gammie (2001|) finds analytically or ~ ^ » 0.022 and nu- 
merically a X 0.0247 (see his Sect. 4.1). By using a = 0.01 we un- 
derestimate his critical accretion rate for fragmentation by a factor 
of order unity. Given the fact that the true critical Toomre parame- 
ter is still under debate (cf. Sect. l5.lt . using a = 0.01 seems to us 
to be adequate. Using a = 0.3 for marginal selfgravity, as Tan & 
Blackman do, rises the critical accretion rate by a about one and a 
half orders of magnitude. 

The role of selfgravity in Pop iii is different compared to those 
found in Pop i. Usual low-temperature dust dominated Pop i discs 
become vertically selfgravitating (KSG) as the disc mass exceeds 
the fraction - of the accreting mass and become fully selfgravi- 
tating (FSG) for disc masses larger than the accreting mass while 
remaining vertical selfgravitating (Duschl et al. 2000). However, in 
the Pop III case the outer parts can be FSG while not any more being 
KSG due to the local mass density. 

This can be understood by examining the ratio 



7.2.2 Thermal and viscous instabilities 

In Fig.|2|we already indicate thermal instabilities. We now want 
to examine the effect of these instabilities by varying the accretion 
rate M(S) at given radial distance as shown in Fig. |4] According 
to eq. J13> the disc is viscously unstable for negative M-I, slope. 
This is the case at very high accretion rates (Lightman-Eardley in- 
stability) in both chemical compositions and also at some kinks in 
the molecular opacities of Pop i material. In our Pop iii disc models 
we do not find any viscous instabilities except at the turning points 
between the thermally stable and unstable solutions. 

Fig- ID shows that limit cycles are present. For both POP i and 
POP III matter the limit cycles are mainly due to H" -absorption. 
The influence of molecular lines (TiO, CO, etc.), however, causes 
further kinks at lower accretion rates and prevents Pop i disks from 
undergoing large AM - AS outbursts. For POP iii no such molecu- 
lar lines are present at these temperatures (H2 lines set in at lower 
temperatures). At the transition from the intermediate CIA cooled 
region of the disc to the isothermal outer disc there is an additional 
thermal instability present. This instability arises because at this 
temperatures H2 dissociates while the contribution of atomic hy- 
drogen to the opacity is not yet that large. 



gKSG 



4npr^ 

M 



For the Pop i case, due to cooling we have the decrease of the disc 
H/R ratio (H/R oc R^^, ^ x 2). With assuming the surface density 
S cc R~(, ^ !s 0), the density is p ex 1,/H tx R(-<-'. The ratio 
8ksg/8c.z "^ R^^i*- increases with radius provided ^ - ^ > -2 which 
is mostly satisfied in the Pop i case. 

The scaling in Pop iii is different. H/R oc R and "L oc Ri (see 
Sect. l7.1.H . Hence p ex "L/H oc Ri. Once the disc reaches the 
temperature-locked region (Sect. lTTOl . gaso/gcH '^ f^^ ■ The ver- 
tical selfgravity is diminishing, while the disc mass in the isother- 
mal part still increases. (Mjisotheimai = jlnXrclr oc rJ). 

Thus in the presence of Pop iii matter one needs to distinguish 
between the local KSG-criterion (sksg/ScM ^ 1) and the global 
FSG-criterion M,, ^ M. 



7.3 Timescales 

For the assumptions in our disc model (hydrostatic, thermal and 
chemical equilibrium, 1 + lD one-zone approximation) we need 

Tvisc » T,h » Tz » Th, 

to be satisfied. We show the timescales for sample models in Fig. 
|5| The chemical equilibrium is implicitly assumed to be fullfilled 
everywhere except where we have a significant abundance of H2. 
In this domain we plot the H2 formation timescale. Otherwise there 
are only ionic and/or ion-neutral chemical rections which proceed 
on timescales much smaller than timescales present in the disc. 

The assumption of chemical equilibrium is satisfied even in 
the Pop m case. The H2 formation timescale is orders of magnitudes 
smaller than the vertical timescale except in the outer regions of the 
disc. As long as the density is high enough, the 3-body reactions 
rapidly convert atomic into molecular hydrogen. In the outer parts 
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Figure 4. Accretion rate versus surface density at radial distance of r = 0. 1 AU in disc models of Pop m (left) and Pop i (riglit) material witli central mass of 
IMq. Tlie meaning of tlie lines, dots, etc. is the same as Fig. 121 




i 1 2 3 2 I 

Figure 5. Timescales for Pop iii disc models with accretion rates 10"'' (left) and 10"* (right) Moyr"'. 



of the disc, where th, increases, the strong density dependence of 
the 3-body reaction time scale can be seen. 

However, our H2 formation timescale should only be used 
with care. Our assumption is that this timescale equals the chemical 
reaction timescale. In general we only could assume that Th, is a 
lower limit to the chemical reaction timescale as the consumed H 
still has to be depleted to reach chemical equilibrium. 

As, however, th, is several orders of magnitude smaller than 
the hydrostativ timescale, atomic hydrogen has got time enough to 
be sufficiently depleted not influencing the opacity and therefore 
the disc structure. In the low-temperature region the only important 
absorption mechanism influencing the opacity is CIA. 

The global maximum of the absorption coefficient normalized 
to standard number densities approximately coincide at the same 
level and the same frequency. Due to the quadratic dependence of 
the CIA on the density of the contributing species, once H2 has been 
formed, CIA of H2/H2 pairs is the dominant absorption mechanism 
at high densities. 



8 DISCUSSION 

We have shown the impact of the different chemical composition 
on stationary accretion discs, comparing Pop iii and Pop i material. 
This is mainly due to the lack of dust and heavier molecules (H2O, 
TiO, etc.) absorption at low temperatures. 

While the inner (hot) parts of Pop iii discs do not largely dif- 
fer from Pop I discs due to the dominant influence of hydrogen 
and helium species on the absorption. Pop iii are optically thin and 
isothermal in their outer parts because of temperature locking due 
to H" -absorption in the absence of dust. This leads to the break- 
down of the thin-disc approximation as H/R increases proportional 
to r J . Pop I discs do not have this locking mechanism. Therefore 
H/R decreases with the onset of the vertical selfgravity. This con- 
siderable flaring may lead to a disc structure modified by irradia- 
tion from the accreting star and the inner disc regions. Such effects, 
however, are beyond the scope of this paper. 

The disk models considered in this paper do not take into ac- 
count effects due to the inner or outer boundary condition. Thus 
for protostellar evolution the innermost solutions will be nonexis- 
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tent being part of the protostar, while the outer disk regions may be 
hidden in spherical accretion still ongoing (see TM04) . 

At the very outermost parts of Pop iii accretion discs the ther- 
mal timescale becomes larger than the viscous timescale as H/R 
reaches unity (see Fig.|5}- However it remains to be shown wether 
the breakdown of the thin-disk approximation (Sect. lT.l"!?! remains 
valid when taking into account irradiation. 

The assumption of chemical equilibrium is well satisfied for 
all regions with reasonable disc solutions in contrast to primordial 
protostellar collapse models. Given the much smaller density, the 
H2 formation timescale exceeds the dynmical time scale there. 

Pop III discs have intermediate regions, where they are able 
to cool below the 3440 K of the outer temperature-locked zones, 
where CIA is the dominant absorption effect. The minimal tem- 
perature achievable decreases with increasing accretion rate. For 
accretion rates high enough these discs do have a vertically self- 
gravitating ring (see Fig.|5J- 

This selfgravitating ring may be subject to change after inclu- 
sion of line absorption in the opacity calculations. At low density 
and low temperature we expect primarily H2 quadrupole line ab- 
sorption which may prevent the disc from reaching the isothermal 
outer part. As the molecular lines are narrow and weak at low den- 
sities, the integrated absorption coefficient added to the continuum 
Planck opacities will become constant with decreasing density for 
given temperature(/<:p(p, T) oc T'^), while it does not in the con- 
tinuum case (Kp(p, T) oc pT^). In this line absorption dominated 
regime, the Pop iii opacities may exhibit the same behaviour as Pop 
I opacities do in their dust and ice dominated regime, although at 
much lower absolute values. 

In this case, there may be at least two coexisting solutions: 
The isothermal part will remain (temperature is only fixed due to 
H^ absorption which will remain the dominant absorption in this 
temperature range) with the addition of a stable low-temperature 
solution branch. 

As mentioned above, some authors conclude that Pop iii 
stars must have been very massive. There have been many pro- 
posals of feedback p henom ena which may halt the accretion 
jOmukai & Palla 20 0lLl2003h . Our stability examinations show a 
direct limiting process. While with Pop i discs the accretion rate 
due to the thermal instability can be increased by 1.5 orders of 
magnitude in contrast to nearly 3 orders of magnitude for Pop iii 
discs. Pop III discs are likely to exceed their Eddington Limit quite 
often and therefore produce massive outflows. However it remains 
to be shown in time-dependent simulations of Pop in discs how this 
instability acts. 

The Toomre parameter of the isothermal outer region becomes 
lower than unity if the accretion rate exceeds 1.07 • lO^^Moyr"' 
(S ect.rm. While Tan & Blackman l20ol) use a = 0.3 (foUow- 
ing lGoodmaaj2003h) we use a = 0.0 1 as the value for the viscosity 
parameter following lOammid J2OOII) . 

Discs well above this threshold accretion rate are certainly 
prone to fragmentation. Discs with lower rates can have Toomre 
parameter values 2 < 1 in the intermediate regime but Q > \ in 
the isothermal outer parts. This does not lead to framentation, as 
the viscous timescale in the outer parts of Pop iii accretion disc de- 
creases until the onset of the isothermal part. This means that the 
first stars can grow very large and that it is exceedingly hard to form 
Pop III planets. 
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APPENDIX A: THE IMPLICIT EQUATION FOR THE 
VISCOSITY 

In the case of the o-- viscosity the viscosity can be written as a func- 
tion of density p and temperature T. We have 

P 



dp 



kT AcrT* Ik^M 4 



Q. ^mp 3cpTgff(v) \ 67TV 3) 
making use of the angular momentum equation and 



Tef£ 



with 



krM 4 4nv 
drtv 3 KfM 



6n 



and further setting 



kT 
fimp 



4 4;r 

T2= - Ts = r 

3 KpM 






we get 



Cg + Cr 



Ti + T2V 



- Ti + T2V + TiV^ 



Tiv^ + {T2 - CoCgTi) v^ + (Ti - C0T2 (Cg + Cr)) v 

-Cori(Cg + Cr) = 
This equation has to be solved for given density and tempera- 



ture. 



APPENDIX B: THERMAL AND VISCOUS STABILITY OF 
RADIATIVELY SUPPORTED DISCS 

In this section, variables like x denote logarithmic variables. 

We start with the equations for energy and hydrostatic equi- 
librium (see. eqns.|5|and|3 

9 



-vnn' 



kT 4(r 



P + ^ r rR + 



4o- 
3cTeff 



G\n + 



M 



jjnip 3cTgjf \ "^ 3 / \ 4pr^ ) 

Using eqns. J4j, Js} and J5J we can write p and T as a function 
of Z. V and r^g remain functions of p, T and S. 

p(E) r(I) Teff(p,r,E) v(p,T,l) (Bl) 

Bl Viscous instability criterion 

We are interested in the derivative 



dM 



dT 



dM dp dMdT dM ^^^^ 

^ + ^ ^ + ^ (B2) 

dp dl. dT dl. dl. 



This derivative has to be calculated for equilibrium pressure 
and energy. Thus we have 



d{Q'-Q-) = 
d[Pg + Pr-g,,nZ + gKS.GZ) = 



which results in 2 linear equations for the unknown derivatives -^ 



a-p 



andf 








.dp dT ^ 
A^ + B^ + C 
dl. dl. 


= 






dp dT 

D-& +E^ +F 
dl dl 


= 






which can be calculated for AE - 


-BDi=Olo 




dp 


BF-CE 


dT CD- 


-AF 


dl 


AE-BD 


dl AE- 


-BD 



With a further examination of eq. <B2> we get 



dM 



dl 



dM dp dM dT dM 
'dpdt Wdl "dl 
dv dp dv dT dv 

a^af '^df'di ^ af "^ 

f[A,B,C,D,E,E,^) 
AE-BD 



dM 

af" 



f[A,B,C,D,E,F,^) 



AE-BD 



<0 



B2 Thermal instability criterion 

Thermal instability requires 



dff_ 
dT 



dQ_ 
dT 



>0 



to be satisfied. We expand Q^ and Q in terms of p and T 



d{Q*-Q-) 



dQ* dQ- 
dp dp 



dp + 



dQ^ dQ- 
dT dT 



dT 



(B3) 



(B4) 



= Adp + BdT 

We further need to assure dP{p, T) = 

d{P,.H+PKSG-Pg+Pi) = Ddp + EdT = 

This is a condition which expresses dp in terms of dT. We get our 
thermal instability condition according to eq. <B4> to 



dff_ 
dT 



dQ_ 

dT 



AE 

-i-fi>0 

D 



AsD 



dp 



f >0 



AE-BD<0 



(B5) 



The term AE - BD is the same as in the denominator of the vis- 
cous instability criterion. Therefore it can be conjectured that if the 
disc is viscously unstable there also may be a thermal instability 
present. This leads to the expression of thermal-viscous instabil- 
ity. However this relation does not hold in a strict sense. There are 
at least formally possibilities having viscous instabilities without 
having thermal instabilities. 

The conditions <B5> and <B3> can in principle be expressed 
now depending only on Ar/p, Br/p, fip and i/ksg with contributions 
of the p-gradients in density and temperature. However we restrict 



ourselves only to supply handy formulae for relevant limit cases. 
(See SectionlOandl?^! 



APPENDIX C: ESTIMATION OF THE H, FORMATION 
TIME SCALE 

We calculate the disc properties analytically for a selfgravitating 
ring, assuming fSp = I. 

^From the angular momentum equation Q with the our vis- 
cosity prescription JF) we get 

kT I r3 
^irtp \ GM 

This leads to 



X(r): 



SnakT 



Plugging this into the equation for the hydrostatic equilibrium leads 
to0 

kT _ I M \ M^fi-ml GM 
^^ " ^ i V^ "^ ""j 9n'-a^k^T^ ~' 

This is an equation which can be solved for T(p). 

+ n\ 



M 



Apr^ j 9n-a-k^p r^ 

2010 K ■£{M,p,r)-(-^\ ^ , , 

^ ^ \0.0l/ I in-4M_„r-i 



lO-'Moyr- 

5 / ,- \-l /,,\1 



10"^ g cm 



-3 



(^)(?) 
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APPENDIX D: COMPARISON TO THE TAN & MCKEE 
(2004) DISC MODELS 

ITM04| present disc models for the inner regions of the accretion 
disc around a Pop in star. Their disc models span 1 decade in ra- 
dial distance (with the exception of the 100 Mq case), are strongly 
influenced by the inner boundary condition. In order to compare to 
our model, we subjected our discs to the same inner boundary. We 
present our results for two of their parameter sets in Fie. lCll 

The results are very similar at the high-temperature region. 
However, the TM04 neglect the low-temperature region. Our mod- 
els show that there are two stable solutions coexisting over large 
radial distances, connected through the thermally unstable part (cf. 
Fig. 121, giving rise to thermal instabilities. 

|TM04 further mention the importance of the ionisation term 
in the energy equation, while the thermal term is negligible. We do 
not account for this effect, as we want to keep the presentation as 
general as possible (see Sect. 17. 1.11 . 

If there is any important term in the energy equation other 
than viscous dissipation and radiative losses, the discs are no longer 
geometrically thin disks. The hydrostatic equilibrium ^ reads 

fj2 Ql A A 



f- a:- + AnGp 



■ 7 • — = (1 - '7ksg) • 7 • 



where u^ = f2r is the azimuthal velocity, y the adiabatic index, and 
P = ypc\ the speed of sound. Disks are thin (HjR « I) as long 
as they are cool (thermal energy small compared to kinetic energy 
of the azimuthal flow). Disks become hotter to reach the same HjR 
ratio either for 7 < f (ionisation or dissociation present) or for 
J/KSG > (see. HjRir) in Fig.|2|for the isothermal outer part of 
the lO^^'Moyr"' and lO^^Mgyr"'). The first effect is present in the 
TM04 models whereas we set y = const. However TM04 do not 
include vertical selfgravity which is the dominant vertical force at 
these high accretion rates. 



with 



f(M,p,r)= 1-1-4.76 



M 



Mo/\10-'*gcm-3 



iM 



Whhk4 = 1.83 • 10-3' (ji^) cm^ s"' JPallaet alllQsl andp = 
/iinpriH we get for the H2 formation time scale: 



Th, 



1 



2k4njj 

7.65 ■ 10^- s • (—^] 
V300K/ 



P 



10^* g cm~3 

M 



ri) 



0.51s. ^(M,p,.).(^f(^;^^^^_,^ 

Mo/ ^10-8gcm-3J I AU/ Vl/ 



Note : For f = 1 we are in the vertical selfgravitating domain, f > 1 
elsewhere. 

Alternatively we could have taken the reaction rate ky, = 
1.3 • 10-32 (^;j;^) ' cm' S-' of lAbel et alj J2002h . Then the cor- 
responding numerical factors would be 1.07 and 7.23. Using k^b 
instead of k4 leads to an order of magnitude higher value of Thj . 

However, with more and more H2 being produced, another 3- 
body reaction (k(, of Palla et al. 1983*) will help in further reducing 
the formation time scale. 
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Figure CI. Accretion disc models for two cases of Tan & McKee (2004): M = IOMq, r, = 300Rq, M = 6.4 ■ lO^^Moyr"' (left), M = IOOMq, r 
M = 2.4 ■ 10^'Moyr"' (right), a was chosen 0.01. We discard all solutions where S.s^, an approximation for the disc mass exceeds the central mass. 
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